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Abstract 

To obtain improved parton densities of the proton, we present a new global analysis 
of deep inelastic and related data including, in particular, the recent measurements of F2 
at HERA, of the asymmetry of the rapidity distributions of production at the FNAL 
pp collider and of the asymmetry in Drell-Yan production in pp and pn collisions. We also 
incorporate data to determine the flavour dependence of the quark sea distributions. We 
find that the behaviour of the partons at small x is consistent with the precocious onset 
of BFKL leading log(l/x) dynamics. We discuss the ambiguities remaining in the gluon 
distribution. We present improved predictions for W boson (and t quark) production at 
the FNAL pp collider. 



1. Introduction 

The increased precision in the experimental measurements of deep inelastic scattering and 
related processes over the last few years has led to a considerable improvement in our knowl- 
edge of the parton distributions of the proton. However several significant sets of measurements 
have become available since the last global analyses of the data were performed to determine 
the parton densities. The new data may be divided into two groups. First we have the mea- 
surements of the structure function F2{x,Q'^) for electron-proton deep-inelastic scattering for 
X < 10^^ by the HI and ZEUS collaborations at HERA ||, 0], and secondly the measurement 
of the asymmetry in Drell-Yan production in pp and pn collisions by the NA51 collaboration at 
CERN 0], and of the asymmetry of the rapidity distributions by the CDF collaboration 
at FNAL 1^. The two groups of measurements have quite distinct implications for the partons 
and are, therefore to a great extent, decoupled. The HERA F2 data offer the first constraints 
on partons in the previously unexplored small x regime, whereas the two asymmetry measure- 
ments probe fine details of the quark distributions in the region x ~ 0.1. The latter information 
is crucial, for example, for a precise determination of the mass of the W boson at the FNAL 
collider. 

Fig. 1 shows the dramatic rise of F2 with decreasing x, which was first observed by the HI 
and ZEUS collaborations [|I|. Also shown are the extrapolated predictions for F2 at = 15 
GeV^ from two parton sets (Dq,D'_) |^ made before the HERA measurements became available, 
but which it was believed would span the data. The upper limit, the curve D'_, was motivated 
by assuming the precocious onset of BFKL dynamics in which the gluon and sea quark 
distributions have the singular form 

xg, xgsea ~ (1) 

as X — i> 0, whereas the Dq curve was the lower limit anticipated from conventional Regge 
expectations with 

xg, xgsea ~ (2) 

Such extrapolations are notoriously unreliable and have failed in the past. Moreover, as we 
will see in section 3, the connection with precocious BFKL behaviour turns out to be more 
subtle than (1) would suggest. Although Dq and D'_ are the most recent published MRS sets 
of partons, a set MRS(H) was subsequently made available,^ which was obtained from a global 
analysis that incorporated the measurements [Q] of F2 from the 1992 HERA run. The curve 
denoted by H on Fig. 1 is an example of the quality of the fit. The main feature is that the 
HERA measurements required a small x behaviour of the form 

xg, xgsea ~ (3) 

Higher statistics measurements of F2, obtained from the 1993 HERA run, have just become 
available in preliminary form |^. We incorporate these data in our new global analysis and 
in section 3 we discuss the implications for the small x behaviour of the partons and for QCD 
dynamics. The new HERA data for F2 are in line with the old (i.e. show the same rise with 
decreasing x), and in fact the MRS(H) partons still give an excellent fit in the HERA small x 
region. 

brief description of MRS(H) can be found in Ref. pO[ |. 
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The data in Figs. 2 and 3 are respectively the measurement of the asymmetry in Drell-Yan 
production in pp and pn coUisions 



Ady = (4) 

and of the asymmetry of the rapidity distributions of the charged leptons from — > £^z/ 
decays at FNAL §, 

- aji-) 

^^^^^ = a{n + a{i-y 

In (1), a = (P a / dM dy\y=Q where M and y are the invariant mass and rapidity of the produced 
lepton-pair, while in cr(^^) = da/dyp are the differential pp —>■ W'^X —>■ l^vX cross sections 
for producing leptons of rapidity y^. Also shown are the predictions of MRS(H) and the 
equivalent set of partons, CTEQ2M, obtained by the CTEQ collaboration |]TT|. We see that 
neither set of partons gives a satisfactory description of both asymmetries. This deficiency 
of the parton sets is not surprising. The reason is that the high-precision muon and neutrino 
deep-inelastic structure function data, which provide the core constraints of the global analyses, 
do not pin down the combination d — u oi parton densities. Indeed the Drell-Yan asymmetry 
experiment was proposed [|12| as it was uniquely equipped to determine just this combination of 
densities. The asymmetry data therefore offer a fine-tuning of the u, d, u and d parton densities 
in the region x ~ 0.1, which is invaluable for the precision studies of the W boson at FNAL. To 
this end we include for the first time the asymmetry data in the global analysis (together with 
the new HERA measurements of F2) and find a new set of partons, which we denote MRS (A). 
The resulting description of the asymmetry measurements is also shown in Figs. 2 and 3. We 
discuss this aspect of the global analysis in section 4. 

The outline of the paper is as follows. We first explain, in section 2, the procedure that we 
follow to determine the parton densities from a global analysis of the data. We give details of 
the new improved MRS (A) parton distributions, and we compare them with the MRS(H) set. 
Sections 3 and 4 consider the impact of the new small x and asymmetry data respectively. In 
section 5 we discuss the ambiguities in the present knowledge of the gluon. In section 6 we 
update the predictions for W boson and top quark production at the FNAL pp collider, and 
finally in section 7 we present our conclusions. 



2. The global analysis 

The parton distributions fi are determined from a global fit to a wide range of deep-inelastic 
and related data. The basic procedure is to parametrize the fi at a sufficiently large Qq {Q1 = 4 
GeV^) so that fi{x,Q'^) can be calculated at higher in perturbative QCD using next-to- 
leading order Altarelli-Parisi (GLAP) evolution equations. In view of the quantity and variety 
of data that are fitted, it is remarkable that an excellent description can be obtained with the 
following simple parametrization 

= AuX^^ (1 - x)"^^ (1 + e„yx + 7„a;) 

xd^ = Aax'^-'il-xy^il + eaVx + ldx) 

xS = Asx-\l-x)'''{l + esVx + lsx) 

xg = Agx~^{l-x)''^il + ^gx), (6) 
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where the valence distributions = u — u and = d — d, and where the total sea distribution 
5" = 2{u + d + s + c) . We assume that s = s. At present there are not enough experimental 
constraints on the gluon to justify the introduction of an extra parameter eg in xg, or to 
determine the exponent A independent of that of the sea-quark distribution S. Three of the 
four Ai coefficients are determined by the momentum and flavour sum rules. The distributions 
are deflned in the MS renormalization and factorization scheme and the QCD scale parameter 
Ajj^{nf = 4) is taken clS 8b free parameter. 

The flavour structure of the quark sea is taken to be 



2u = 


0.4(1 


-6)S-A 


2d = 


0.4(1 


-6)S + A 


2s = 


0.2(1 


-6)S 


2c = 


6S 





(7) 

at = Ql = A GeV^ with 

xA = x{d-u) = A^x'^'il- xy%l + -fAx). (8) 

The flrst hint that the u, d flavour symmetry of the sea is broken (with ci > m on average) came 
from the evaluation of the Gottfried sum by NMC []T3[. In fact until then all global analyses 



had assumed u = d. As we shall see in section 4, the observed Drell-Yan asymmetry provides 
further evidence that d > u, which we allow through the parametrization of A given in (H). 
The flavour breaking can be associated with the breaking of p — 02 meson Regge exchange 
degeneracy and so we choose the exponent rji in (|^) to be the same as that of one of the valence 
quark densities. 

As in earlier MRS parametrizations, we assume that the input strange sea is suppressed 
by 50% in relation to the u and d sea distributions - hence the factors 0.4, 0.4 and 0.2 in 
Eq. (0). The first indications for such a suppression came from early deep inelastic dimuon 
production data, but now there is much firmer evidence for our input assumption. The CCFR 



collaboration WM have performed a next-to-leading order analysis of their uN /i~/i+X data 
and deduce that the strange sea distributions should lie within the shaded band shown in Fig. 4. 
The strange sea that we find is shown by the curve denoted by MRS (A) in Fig. 4 and satisfies 
the experimental constraint. 

The input charm sea is determined by the EMC deep inelastic data for the structure function 
We proceed as follows. We assume that 

c(x, = for g2 < ml (9) 

and generate a non-zero distribution at higher by massless GLAP evolution at next-to- 
leading order. The structure functions are also calculated using massless coefficient functions. 
Since = ml falls below our input scale, Ql = A GeV^, we use an approximate set of partons 
to evolve between and 4 GeV^. Taking zero charm at = we find that the shape 
of the resulting charm distribution generated at = 4 GeV^ is well described by the input 
parametrization of the overall S distribution. As might be expected the normalization, specified 
by the parameter 5 of (|^), depends sensitively on the value chosen for mc- We adjust 6 to give a 
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good description of the data for Q >'o GeV^. The fit is shown in Fig. 5 and corresponds to 
5 = 0.02, and to a charm distribution which can be evolved from zero at = ml = 2.7 GeV^. 
The value 6 = 0.02 implies that, at the the input scale Qq = 4 GeV^, the charm sea carries 
0.4% of the proton's momentum, as compared to nearly 4% carried by the strange quark sea. 

Our prescription for the charm quark distribution is only valid far above threshold, = 
— x)/x 3> 4m^. Near threshold, a more rigorous treatment of quark mass effects is 



required. As discussed in detail in Ref. (see also ITg), various prescriptions are possible. 
For example, one can define heavy quark densities according to our prescription, and absorb the 
threshold effects into the coefficient functions. Nevertheless, Fig. 5 shows that our treatment 
of the charm distribution does give a reasonable description of the EMC data [|I^ for F2 with 
> h GeV^. These data do in fact lie in the region W"^ ^ 4m^. It is interesting to note 
that our satisfactory description means that these EMC data show no necessity for a small 
"intrinsic" or non-perturbative charm component, as advocated by Brodsky et al. except 



possibly for one data point at x = 0.42 (not shown) which lies well above our fit. Also shown 
in Fig. 5 are the predictions in the small x region which are relevant for future measurements 
at HERA. Note that small x implies large W^, and so our treatment should be reliable in this 
region. 

It has been argued [|19] that the strange sea as measured in neutrino scattering should be 
different from that in muon scattering on account of the different mass thresholds in W*g —>■ sc 
as compared to 'y*g — ss. In practice the neutrino data have been corrected by the CCFR 
collaboration to approximately take into account the rric 7^ effects and to allow for this 
difference. Our strange and charm seas should be interpreted as those for muon data. Our 
treatment of the charm sea at large x and modest is suspect, but then its contribution is 
too small to distort the analysis. Finally, the 6-quark contribution is included by assuming that 

b{x, Q^) = for < ml (10) 

with ml = 30 GeV^ and non-zero contributions at higher generated by the same "massless" 
prescription that was used for c{x,Q^). 

The experimental data that are used to constrain the parton densities are listed in Table 1, 
together with the leading order partonic subprocesses, which helps us to see which features 
of the distributions are constrained by the various data sets. As far as the global analysis is 
concerned, the "core" deep-inelastic data for x > 0.01 are the BCDMS [0 and NMC [|, || 
measurements of and F2^, and the CCFR measurements of F2^ and xF^^ . It is 
useful to inspect the leading order expressions of these structure functions in terms of the 
quark densities 

F^P-Fr = lx{u + u-d-d) (11) 
KF^ + FD = ^x{u + u + d + d+ls) (12) 
pvN^puN ^ ^(^u + u + d + d + 2s) (13) 
\x{Ff + F^^) = xiu-u + d-d), (14) 

where is an isoscalar nuclear target and where, for simplicity, we have neglected the small 
c quark contribution. On a [x, Q^) point-by-point basis, these four observables determine four 
combinations of quark densities, which we may take to he u + u, d + d, u + d and s. The 
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Process/ 
Experiment 


Leading order 
subprocess 


Pcirton determination 


DIS (/xiV fj,X) 

BCDMS, NMC 

DIS {vN fiX) 

CCFR (CDHSW) 


> 

W*q q' 

> 


> 


Four structure functions — > 

u + u 
d + d 

u + d 

s (assumed = s), 
but only / xg{x)dx ~ 0.5 
[u — dis not determined] 


fj,N ccX 

F|, EMC 


7*c c 




c 0.1s at Ql 


uN ix+n X 
CCFR 


W*s c 




s K,l;u (or ij) 


DIS (HERA) 

F^^ (H1,ZEUS) 


Yq^q 




A 

[xq ^ xg x~^, via g —>■ qq) 


PP ~^ 1^ 

WA70 (UA6) 


qg^iq 




q(x ?a 0.4) 


pN n+fi-X 

E605 


qq 7* 




q= ...{l-x)'^s 


pp,pn — > /x+/Lt~X 

AT A CT 1 


uu, dd 7* 
ud, du — > 7* 




(u-d) atx^ 0.18 


PP WX{ZX) 
UA2, CDF, DO 

— > asym 
CDF 


ud^W 




u, d at X1X2S ~ —>■ 
X ^ 0.13 CERN 
X ^ 0.05 FNAL 

slope oi u/d Bit X 0.05 



Table 1: The experimental data used to determine the MRS parton distributions. The last 
column gives an indication of the main type of constraint imposed by a particular set of data. 
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difference d — u is not determined and, moreover, at leading-order it appears that the only 
constraint on the gluon is through the momentum sum rule. 

From (0) and ([121) we see that the strange quark distribution is essentially determined by 
the structure function difference 

xs{x) ~ iF^^'ix) -SF^'^ix). (15) 

At small X ^ 0.08 the value obtained is significantly larger than that found in the dimuon 
analyses that we mentioned above. We believe the dimuon result is more reliable, as the 
determination based on (|T5|) is very sensitive to the relative normalization of two different data 
sets and to the heavy target corrections that have to be applied to the neutrino data. However 
it does mean that in the global analysis this discrepancy will show up in the description of 
and/or F2^ at small x. The quality of the MRS(A) fit to the deep-inelastic data is shown 
in Figs. 6-8. We note that the discrepancy occurs in the description of F2^ at x < 0.05, see 
Fig. 8(a). We allow an overall normalization parameter for the CCFR data. When we fit to all 
data (satisfying our > 5 GeV^, W"^ > 10 GeV^ criteria) we find a normalization factor of 
0.95. If, however, we were to exclude the x < 0.08 CCFR data then the normalization factor 
becomes 0.97 and the description of the remaining neutrino data improves (x^ = 90/132 data 
points as compared to 171/160). 

Heavy nuclear target corrections are apphed to the neutrino data, and deuteron screening 
corrections are made to the small x muon-deuterium data as described in Ref. [^. As in 
the earlier analyses ^ |T0[ the WA70 prompt photon data |^ and the E605 Drell-Yan 



production data [24| are included in the fit. The former constrain the gluon and the latter pin 
down the shape of the sea quark distributions. 

The values of the parameters of the starting distributions, (^), of the new MRS (A) set of 
partons are listed in Table 2. In addition, the value of the QCD scale parameter is found to be 
Ajj^{nf = 4) = 230 MeV, which corresponds to as(M|) = O.II25, as it was for the MRS(H) 
set of partons. The parameter values for this latter set are also listed in Table 2. Figure 9 
shows the MRS (A) parton distributions as a function of x for two different values of Q^. A 
comparison of the "new" MRS (A) and "old" MRS(H) partons is shown when we discuss the fit 
to the Drell-Yan asymmetry measurement in section 4(a), since the differences between the two 
sets arise mainly from introducing this data point into the global analysis. Finally, we show in 
Table 3 how the proton's momentum is shared among the various parton flavours in the new 
MRS (A) set at different values. 



3. Small X behaviour 

The new measurements of F2^ obtained by the ZEUS and HI p| collaborations in the 
low X regime, x < 0.005, are shown in Figs. 10 and 11. These data were included in the global 
analysis and are, in fact, the only constraint on the parameter A in (^ which controls the small 
X behaviour of the sea xS ~ x^^ (and of the gluon xg ~ x^^). The existing set of partons, 
MRS(H) with A = 0.3, is found to give an excellent description of the new data, and so it is not 
surprising that the new parametrization, MRS (A), has the same value of A. The description of 
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MRS (A) 


MRS(H) 






0.775 


0.777 


Glue 


A 


-0.3 


-0.3 






5.3 


5.3 




Ig 


5.2 


5.2 




Til 


0.538 


0.335^ 






3.96 


3.90+ 




e« 


-0.39 


4.40+ 


Valence 




5.13 


8.95+ 






0.330 


0.224 






4.71 


4.65 




Gd 


5.03 


44.3 




Id 


5.56 


13.2 




As 


0.411 


0.386 




Vs 


9.27 


9.01 


Sea 


es 


-1.15 


0.11 




Is 


15.6 


12.6 




Aa 


0.099 


0.055 




7A 


25.0 





Table 2: The numerical values of the starting distributions (6) of the MRS(A) set of partons. 
For comparison we also list the values corresponding to the MRS(H) partongj. Note that Ag is 
fixed by the momentum sum rule, and is therefore not a free parameter. 



the HERA data is shown in Figs. 10 and 11, together with that of the extrapolations obtained 
from the pre-HERA MRS(D[),D'_) parton sets. 



(a) Predictions from perturbative QCD 

The value of A is of great importance for understanding QCD dynamics in the small x 
regime. Therefore before we discuss the uncertainty and the implications of the A = 0.3 
determination, it is useful to briefly review the perturbative QCD expectations for the small x 
behaviour of F2. In the small x regime we encounter large log(l/x) contributions which have 
to be resummed. It is necessary to distinguish two different limits. First, the BFKL limit of 
small X and moderate > Qq in which we sum the large log(l/x) terms, but keep the full 

dependence, not just the leading logQ^ contributions. Second, the small x and large Q'^/Ql 
region where the double-leading-logarithm approximation of GLAP evolution is appropriate, 
in which the aslog(l/x)log(Q^/(5o) contributions are resummed. Ideally we seek a formalism 



which embodies both limits. Some progress has been made by Marchesini et al. [25] to find a 
unified evolution equation, but much remains to be done before we can extract the form of the 
small X behaviour of partons. However we explain below why this is not a serious obstacle to 
the extraction of partons from the small x HERA data. 



^For MRS(H) and previous MRS analyses we parametrized x{uv +(iv) by the expression that we use here for 
xu^. Thus parameters marked in Table 2 correspond to x{u^ + d^) and not to xu^. To improve the precision 
at small x we have repeated the MRS(H) analysis and so the MRS(H) parameters listed in the Table are not 
precisely the same as those of Ref. |lO[ . 
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(a) I^UreV ) 


9 




rl 

Civ 




LiX 


9e 




9/i 


7 




97 9 


Q 8 


u.u 


o.u 


o.o 


u.o 




20 


44.8 


25.1 


9.1 


6.3 


8.6 


4.2 


1.4 




100 


46.0 


22.9 


8.2 


6.5 


8.6 


4.7 


2.0 


0.7 


10^ 


47.2 


20.5 


7.4 


6.8 


8.7 


5.1 


2.8 


1.5 



Table 3: Fractions (in %) of the total proton momentum carried by the various partons in the 
MRS (A) set. 



First we review the BFKL expectations at small x. The BFKL equation is effectively 
the leading aslog(l/a;) resummation of soft gluon emissions. The equation may be solved 
numerically and the fc-r-factorization theorem used to predict the small x behaviour of ¥2- It 
is found [EBl that 



(16) 



where the coefficient C of the BFKL contribution and the non-BFKL term, F^^ ^ are weakly 
dependent on x. The magnitude of C is dependent on the treatment of the infrared region of 
the integration over the transverse momenta of the emitted gluons - but a physically reasonable 
choice of the infrared parameter yields an ^2(0;, Q^) in good agreement with the HERA data. 
The important point is that the value of the exponent, ~ 0.5, is a stable prediction. That 
is, the x~2 shape, in ([16|) is a characteristic property of (leading-order) BFKL dynamics; is 
not a free parameter but is determined dynamically. (For this reason we use the subscript L to 
distinguish the "Lipatov" A^ from the free parameter A in Eq. (^.) Of course the predicted small 
X behaviour may be modified by sub-leading corrections. At sufficiently small x, shadowing 
corrections will almost certainly suppress the growth of F2 with decreasing x, but this effect 
is expected to be small in the HERA regime, unless "hot-spot" concentrations of gluons occur 



within the proton p6|. The non-BFKL term, Fg^, is one sub-leading contribution which may 
be estimated 



F^'ix 



0.1, 



0.4, 



(17) 



or perhaps rising slowly with decreasing x as implied by the "soft" Pomeron. 

We now turn to the second limit, that is the form of F2 found at small x and large Q'^/Qq 
from the double-leading-logarithm approximation of conventional GLAP evolution. In this case 
the result depends on the choice of the exponent A in the input gluon and sea quark distributions 



X 



-A 



xq^XjQq), xS{x,Qq) 
as X — 0. If non-singular input forms are used with A < 0, then 

F2{x,Q') ~ exp(2[e(gg,Q2)log(l/x)] 
where the "evolution length" 

rQ' dq^ 3«,(g2) 



TT 



(18) 



(19) 



(20) 
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That is F2 increases faster than any power of log(l/x) but slower than any power of x. We see 
that the steepness in x is not stable to evolution in Q^, but increases rapidly with the evolution 
length ^(Qg, Q"^)- On the other hand if singular input forms, (|18|), are chosen with A > 0, then 
the x~'^ shape 

F2{x,Q^) ~ h{Q^)x~'^ (21) 

is stable to evolution in Q^, and the behaviour ( pi]) overrides the double leading logarithmic 
structure of (111) p^. 



(b) Implications of the global analysis in the HERA regime 

Can the HERA measurements of F2 distinguish between the small x behaviours presented 
in (p!6D, ( p!9D and (|2T|), bearing in mind that the data are well described with MRS (A) partons 
with a sea quark input form 

xS{x,Ql) ~ x-^{l + tsV^+ ...) (22) 

with A = 0.3, and that F2 mirrors the behaviour of xS*? 

It is informative to first recall the GRV dynamical parton model of Gliick, Reya and Vogt 



28| , which pre-dated the HERA small x measurements. GRV predicted a steep behaviour 
of F2, based on the double-leading-logarithm form (|19D, by evolving from valence-like input 
distributions at a very low scale, Qo = 0.3 GeV^. Although the GLAP small x forms of (|T|) 
and (|21| ) are quite distinct, in a limited (x, Q"^) region about, say, (x, Q^) the double-leading- 
logarithm form mimics an x~^ behaviour with 

. _ /36 1og[log(gVA^)/log(gg/A^)] \^ 

where, for five flavours, 60 = 23 and A 150 MeV. If we were to take = 0.3 GeV^, as in 
the GRV model ||28|, then in the HERA regime we have A ^ 0.4. As we have seen, the HERA 



data appear to prefer a slightly lower A. We note that the observed x~°'^ behaviour can be 
approximately mimicked if we were to evolve from non-singular input forms at a higher starting 
scale of Qo = 1 or 2 GeV^, see, for example, Ref. [^. 



The detailed x and behaviours of the GLAP-based forms, (IH) and (^T]), are quite 
distinct. However if Qq is taken as a free parameter in ([T9|) and A as a free parameter in (|2l|) 
(or rather in (^2])), then the HERA data are not yet precise enough to distinguish between 
them. It will be even more difficult to distinguish between the BFKL form ([T6|) and the GLAP- 
based behaviour given by (pT]). In principle it might appear that the different dependence 
of C(x, Q"^) ~ (<5^)2 in ([16|) , and of h^Q"^) ~ logQ^ in ([2l|), would be a sufficient discriminator, 
but in practice the differences are not large |3^ . 



The most distinctive theoretical prediction of section (b) is the behaviour with A^ ~ 
0.5, embodied in (|1^), which was obtained from the leading-order BFKL equation. Since A^, is 
not a free parameter, the HERA data could have ruled out the precocious onset of this BFKL 
behaviour. Rather the data are found to be consistent with (0), which may be written in the 
form 

F2 ~ x"5(l + ax5) (24) 
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where a = F^^/C is weakly dependent on x. 

At Q2 ~ 15 GeV^ we estimate! that a ^ 30. 
Suppose we approximate the BFKL expectation, (|2^) , by the simpler form F2 ~ then it 
follows that ^ 

A ^ 1 ~ 0.3 (25) 



2 2(1 



ax 2 



if we insert a 30 and a typical value of x in the HERA regime, x ~ 5 x 10"^. That is, 
the value of A predicted from BFKL dynamics is in agreement with that found in the global 
analysis. This very approximate identification at = 15 GeV^ was simply made to illustrate 
the similarity of the BFKL description with the global fit based on GLAP evolution. In fact 
present data for F2 cannot distinguish between the two descriptions. An indication of the 
accuracy required to discriminate between them can be found by comparing the MRS(H) and 
AKMS curves in Fig. 3 of Ref. g. 



To identify the BFKL behaviour (^) with (p2D we have used the fact that the es term gives 
a small contribution for x ~ 10~^ in the MRS (A) fit. It will be no surprise to report that in 
the global analysis there is a strong correlation between the values of the A and parameters 
in (p2D. Almost as good a global description of the data can be achieved with, for example. 



A = 0.2, es = -3.3 (26) 

or with 

A = 0.4, es = 2.9 (27) 

In other words, the value of the parameter A in Eq. (|^) is not well determined by the HERA 
data. 

It is important to note that it is the sea quark densities, and not the gluon, which are 
constrained by the HERA measurements of F2 at small x. Since the sea quarks are driven by 
the gluon, via g qq, we have assumed that they have a common behaviour at small x, 
see Eq. (P). However there is, as yet, no experimental confirmation of this assumption, and 
the ambiguity in the gluon distribution is by far the largest uncertainty in the parton densities. 
It is therefore crucial to make a direct determination of the gluon in the region x < 0.05. In 
section 5 we discuss the gluon distribution in more detail. 

To summarize, we find that, within the global analysis, the HERA measurements of F2 are 
well fitted by a parametrization which embodies a small x behaviour of the form F2 ~ xS ~ 
^-0.3 'jj^jg behaviour is consistent with the expectations of BFKL dynamics, as expressed by 
F2 ~ + of ([T6|) , since this form mimics an x~^'^ behaviour in the HERA x regime 

under consideration. The observed x~° '^ dependence is however not consistent with the GRV 
"valence" model which gives a steeper behaviour in the HERA region (with an effective A > 0.4). 
Such an approach could be made consistent with the present data if a higher starting scale Qq 
were used to decrease the evolution length. We also noted that very precise measurements of 
F2{x,Q'^) at HERA will be needed to distinguish between the BFKL behaviour and GLAP 
evolution from "singular" input distributions (motivated by BFKL dynamics). The converse 
of this result is of great practical benefit to the determination of partons. It means that we 



^The most recent AKMS Q description, @, of HERA data (which can be found in Ref. [|o|) has the 
approximate form F2 = O.OlSa;"^ + 0.385{x/xo)~°-°^ at = 15 GeV^ where xq = 0.1. Thus we have 
a = F^^/C w 0.6/0.018 w 30 for x in the HERA regime (x - 5 x 10"''). 
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can base the global analysis of the data, including the small x measurements of F2, on GLAP 
evolution. When precise data for F2 become available over a range of x and in the HERA 
regime, a global analysis may show a systematic departure from the GLAP forms which would 
indicate BFKL dynamics and/or the onset of shadowing corrections. 



4. Constraints from the asymmetry measurements 

As we saw in the introduction (Figs. 2 and 3) the existing set of partons MRS(H), and its 
counterpart CTEQ2M, do not simultaneously describe the Drell-Yan and W rapidity asymme- 
try measurements. This is one of the major motivations to redo the global analysis incorporating 
these data, and to present a new set of partons - MRS (A) - which rectifies the deficiency. The 
improvement is shown by the MRS (A) curves in Figs. 2 and 3. 



(a) Drell-Yan asymmetry 

We have seen, from inspection of ([Tl|)-(p^), that the "core" deep-inelastic structure function 



data offer little constraint on d — u. An indication that d ^ u came from the NMC evaluation 
r3 of the Gottfried sum 



/gs = ['^{Fr-FD = I -/sea, (28) 

Jo X 

where the last equality separates the sum into a valence and a sea quark contribution. NMC 
directly measured the integral over the interval 0.004 < a; < 0.8 and found 0.236 ±0.008 (stat.) 
]T3| . If we correct for deuteron screening effects and assume reasonable extrapolations over the 
unmeasured x intervals then the NMC result implies 

= lj\d-u)dx ^ 0.1 (29) 

for ~ 4 to 7 GeV^. This constraint on d — u given by (pQ]), which is sensitive to the 
assumption made for the small-x behaviour of the valence quarks, is not directly included in 
the global analysis. It is already correctly subsumed since we fit to all the component NMC 
data. 

The asymmetry in Drell-Yan production in pp and pn collisions offers a direct determination 
oi d — u [O . The sensitivity can be seen from the leading-order expression 



^ ^ (ypp - (ypn ^ {u-d){u-d) + ^{uu-dd) 

(Tpp + CTpn [u + d){u + d) + \{uu — dd) + |(ss + 4cc) 



where a = d'^a/dMdy\y=Q and where the partons are to be evaluated at x = ^/T = M/y/s] 
M and y are the invariant mass and rapidity of the produced lepton pair. Since u > d in the 
proton, the asymmetry is positive for parton sets with d — u zero or small, but becomes negative 
a.s d — u increases. 

When we include the recent NA51 asymmetry measurement 0, Ady = —0.09 ±0.02 ±0.025 
at X = 0.18, in the next-to- leading order analysis, the description of the data is shown by the 
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MRS(A) curve in Fig. 2. Fig. 12 shows that the introduction of the NA51 measurement leads to 
a significant increase in J — -u in going from the "old" MRS(H) to the "new" MRS (A) partons. 
In the region 0.02 < x < 0.7, where a full set of deep-inelastic data exist, we expect the MRS(A) 
partons to be very similar to those of MRS(H), except that d — u would be much larger while 
at the same time approximately conserving u + d, d + d and u + u (as expected from (|lT])-(0)). 
Thus we anticipate an increase of d which is compensated by a corresponding decrease in u 
and d, which in turn requires a similar increase in u. The picture is a bit too simplistic since 
there are new NMC data in the MRS (A) fit and moreover we have to maintain the description 
of the W rapidity asymmetry measurements. Nevertheless the comparison of the MRS (A) and 
MRS(H) shown in Fig. 13 displays the trends that we expect. 

The NA51 collaboration also measure the asymmetry, A^, of J/ip production in pp and pn 
collisions. Now J/ ip production can be of gg, as well as qq, origin. We may use a simple model 



31| , in which 

^{gg ^ ^) = r&{qq ^) (31) 

where r = 0.5, to estimate the asymmetry. This value of r is obtained by comparing the rate 
of J/ ip production in pp and pp collisions . To leading order we have 

^ ^ ^pp(^) - ^pn(V') ^ {u-d){u-d) 

app{ip) + a.pn{^) rgg + {u + d){u + d) + 2ss + 2cc' 

where a = da/dy\y=Q and where the partons densities are evaluated at x = M^j ^ and 
= M|. In principle it appears that a measurement of A^ offers a more direct determination 
of (m — d) than Ady, but in practice the gg term dominates and considerably dilutes the 
predicted asymmetry A^. As a consequence, the NA51 measurement 0] = —0.03 ± 
0.002 (stat.) ± 0.02 (sys.) is not able to distinguish between the different sets of partons. 
However the value A^ = —0.013 predicted using the MRS(A) partons is consistent with the 
measured value. 



(b) W charge asymmetry 

The charge asymmetry at the FNAL pp collider 

da{W+)/dy-daiW-)/dy 
^^^^ " da{W+)/dy + da{W-)/dy ^'^'^^ 

is a sensitive probe of the difference between u and d quarks in the x ~ 0.1, Q ~ My/ region. 
Because the u quarks carry more momentum on average than the d quarks, the tend to 
follow the direction of the incoming proton and the W~ that of the antiproton, i.e. Ay/ > for 
y > 0. Thus a precise measurement of the W asymmetry - in practice the asymmetry Ag of 
the charged lepton from W decay (^) - serves as a valuable independent check on the u- and 
(i-quark distributions. In an earlier paper |3^, we showed that there was a direct correlation 
between the lepton asymmetry and the slope of the d/u ratio. To see this, we first note that the 
dominant contribution to production comes from the ud (du) annihilation process. 

Thus 

_ u{xi)d{x2) - d{xi)u{x2) 
u{xi)d{x2) + d{xi)u{x2) 
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where the scale Q = Mw is imphcit for the parton distributions, and 



xi,2 = xo exp(±y) , xo = . (35) 
If we introduce the ratio Rdu{x) = d{x)/u{x), then for small y 

Aw{y)^-xoy^!^^. (36) 
In reality, the situation is of course more complicated - it is the lepton asymmetry (M) which is 



measured and there are subleading and higher-order corrections to Eq. (|34D . Nevertheless, the 
correlation implied by Eq. (|36| ) is evident in the full prediction.! Figure 14 shows the MRS (A), 
MRS(H) and CTEQ2M u, d and d/u parton distributions as a function of x aX = M^. The 
X range is chosen to correspond to the lepton asymmetry measurement by CDF at the FNAL 
pp collider The slope of d/u is significantly larger in magnitude for the CTEQ set, and 
this leads to the larger lepton asymmetry shown in Fig. 3. MRS(H) partons give an excellent 
description of the W asymmetry, in fact slightly better (in terms of x^) than MRS (A) which 
includes the data in the fit. For this reason u and d do not change as much as might be expected 
to compensate for the increase in — u in going from MRS(H) to MRS (A), see Fig. 13 and the 
discussion following Eqs. ([TT| ) - ([T^) . 



From Fig. 3, we can conclude that CTEQ2M is ruled out by the data. This illustrates the 
discriminating power of the CDF data, since CTEQ2M is consistent with all the deep inelastic 
scattering data, and with the NA51 Drell-Yan asymmetry. It is interesting to note that the 
CTEQ2M gives a reasonable description of the NMC F^/Ff data pT| shown in Fig. 7, a quantity 



which is also correlated with the d/u ratio. The reason that there is no contradiction between 
this and the misfit of the CDF asymmetry data is that the n/p ratio is more influenced by the 
u and d distributions than is aw ^ the former is, in a sense, linear in the small q/ q correction, 
^2 ~ |('^* + d) + ^{u + d), whereas the latter is quadratic: a{W~^) ud + du. Hence the W 
asymmetry is a more direct probe of the d/u ratio than the n/p structure function ratio. The 
"incorrect" d/u ratio of CTEQ2M is compensated by the antiquark contributions in the fit to 
the n/p ratio. 

Another place where the d/u ratio is important is in the precision W mass determination in 
pp collisions. Because of the finite rapidity acceptance of the experiments, the lepton transverse 
momentum spectrum - from which the mass is determined - depends on the shape of the 



rapidity distribution of the W p4|. Once again the correlation is with the d/u ratio, and hence 
the W asymmetry can be used to rule out sets of parton distributions in the analysis. In fact 
it should be possible to produce sets of distributions which give "±lcr" fits to the W asymmetry 
data of Fig. 3, and to use these to infer a ±lcr error on Mw from parton distributions.!! We 
shall return to this issue in a future study. 



5. The gluon distribution 

^The curves in Fig. 3 are calculated using the next-to-leading order program DYRAD of Rcf. |Q. We thank 
Nigel Glover for helping with these calculations of the W asymmetry. 

^We are grateful to Arie Bodek and David Saltzberg for discussions on this point. 
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The new data which we have introduced into the global analysis has enabled us to fine-tune 
the quark distributions. In this section, we comment briefly on the status of and prospects for 
information on the gluon distribution. 

The only precise information on the gluon distribution is that the total momentum fraction 
is about 43% at = (Table 3). Traditionally, information on the gluon shape has come 
directly from large px prompt photon production in fixed-target pN collisions, where the qg — *■ 
7g process is dominant. In particular, past and present MRS distributions have used the high 
precision WA70 data [^] to constrain the shape of the gluon in the x ~ 0.3 — 0.4 region. 
Indirect information also comes from scaling violations in fixed-target deep inelastic scattering, 
where the gluon gives an important contribution to the dependence particularly at small 



X. In previous studies |^^, we have investigated the interplay between the prompt photon 
and deep inelastic data, and derived "±lcr" gluons which attempted to span the allowed range 
of the distribution at medium-to-large x. We should stress that all our recent distributions, 
for example MRS(Dq, D'__, H, A), are global "best fits" to the data (available at the time of 
the analysis) of all the processes listed in Table 1. The HERA measurements of F2 were not 
available for the MRS(D) analyses, but these data do not pin down the gluon and so Dq and 
D'_ can be used to demonstrate the ambiguity in its behaviour. However the point we wish 
to emphasize is that the recent global analyses are "best fits" to the WA70 data, and so the 
spread of the gluons obtained underestimates the uncertainty. Nevertheless, it is interesting 
to compare the variation in the recent MRS gluon distributions. Figure 15 shows the Dq, D'_ 
and A gluonsi at = 20 GeV^. We see immediately that there are regions in x where these 
are very different (note the logarithmic scale), and regions where they are quite similar. We 
discuss these in turn below. 



At very smallx, x < 0(0.01), the size and shape of the gluon is controlled by the x~^ term 
in the starting parametrization. Since the gluons are constrained to be similar at larger 
X (see below), a larger A implies a steeper, larger gluon in the very small x region. As 
discussed in detail in section 3, perturbative QCD (whether manifest in GLAP or BFKL 
dynamics) always ties the small x behaviour of the sea quarks and gluons together, and in 
the MRS fits we always use the same A parameter for both, see Eq. (||). Thus the spread 
of the Dq, D'_ and A gluons in Fig. 15 is very similar to that in the corresponding F2^s 
at small x. Note, however, that while we would expect the ratio q/g to be approximately 
constant as x — > 0, we cannot say at which value of x this constancy should set in. In 
fact, for MRS (A) this proportionality between the distributions is valid over a broad 
range of small x, since the term turns out to be comparatively small in the sea-quark 
distribution and since there is no such parameter in the gluon distribution^, i.e. = 
(Eq. i)). 

Now that the very small x quarks are being pinned down by the HERA data, it is impor- 
tant to measure the gluon at small x to look for evidence of a steeply rising distribution. 
Various methods have been suggested for directly extracting the gluon distribution at 
HERA. In the long term, the most promising appears to be the longitudinal structure 



^The A and H gluons are essentially identical, see Fig. 13. 

^Our approach has always been to use simple, physically-motivated forms for the starting distributions and 
to only introduce extra parameters as and when required by the data. 
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function . In the meantime, some very rough indication can be obtained from the 

derivative dF2/ dlogQ"^ which, hke F^, is dominated by a term proportional to the gluon 
at small x. Figure 16 shows xg, and dF2/ dlogQ"^ at small x for the Dq, D'_ and A 



partons. Also shown are the first data on the F2 slope from HI The size of the 

errorbars prevents any conclusions from being drawn at present. 

The production of large px jet pairs with large rapidity at the FNAL pp collider has also 
been suggested as a way of probing the very small x distribution Preliminary 
indications appear to favour a steeply rising gluon, although a more complete study of 
the experimental systematic errors and of the next-to-leading order QCD corrections is 
required before any definite conclusions can be drawn. 

(ii) At medium-large X , 0(0.2) < a; < 0(0.4), the gluon distributions in Fig. 15 are constrained 
by the prompt photon data and are therefore similar. This is not to say that there is 
no uncertainty in the gluon this region - in fact quite a large spread is allowed by the 
errorbars on the data and by the factorization and scale dependence of the next-to-leading 
order QCD cross section [^. Roughly speaking, this translates into a 0(±1) uncertainty 
on the power rjg of (1 — x). Note that in the MRS parametrization of the gluon, Eq. (P), 
the parameter 7^ compensates the different small x behaviours in, for example, Dq, D'_ 
and A, leading to roughly equal distributions with the same parameter rjg at larger x. 
The larger the value of A, the bigger is the corresponding value of 7^. 

(iii) At very large x, x > 0(0.4), there are essentially no constraints on the gluon distribution. 
Phenomenologically, this region is of little importance for high energy hadron-hadron col- 
liders. It is interesting, however, that the power of (1 — x) determined by the prompt 
photon data at medium-large x, rjg = 5.3 for MRS(A) at Qq = 4 GeV^, is perfectly con- 
sistent with the naive expectation from dimensional counting, rjg = 5, which is supposed 
to be valid in the limit x — 1 at some low scale. 

(iv) The medium-small x region, 0(0.01) < x < 0(0.2) is perhaps the most interesting of 
all. Because of (a) the momentum sum rule and (b) the requirement of similar gluons at 
larger x, and gluon distribution which is much larger at very small x must be smaller in 
this region. This leads to a cross-over point at around x ~ 0.01 for the gluons shown in 
Fig. 15. One direct consequence of this is that, in this region of x, structure functions such 
as F2 evolve more slowly with for sets with "singular" gluons at very small x. This is 
illustrated in Fig. 17, which shows at x = 0.05 as a function of Q, with data from NMC 
and curves corresponding to the MRS(Dq, D'_, A, H) and CTEQ2M partons. The fits 
are, of course, constrained to be in agreement in the Q range where accurate NMC data 
are available. Notice, however, the divergence of the predictions as Q increases. This is 
directly related to the magnitude of the gluons in this region of x and Q, as is evident 
from Fig. 18 which shows the corresponding gluons at the same value of x, as a function 
of Q. The effect of the gluon on the evolution of the quark distributions at medium-small 
X has important consequences for W and Z physics at the FNAL pp collider. This will 
be discussed in more detail in the following section. 

There are, in fact, several hard scattering processes at the FNAL collider which offer, at 
least in principle, the possibility of measuring the gluon distribution in the medium-small 
X range. Data on prompt photon production in the range 10 GeV/c ^ < 100 GeV/c 
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has recently become available [50|. Next-to-leading order predictions using "standard" 
gluon sets give a reasonable description of the data for pj, > 30 GeV/c [0. Below this 
value, the measured cross section rises more rapidly than the predictions, indicating an 
excess of data over theory. However this is precisely the region where (a) the uncertainties 
due to scale dependence and to matching the theoretical and experimental definitions of 
"isolated" photons are greatest, and (b) one expects significant contributions from higher- 
order processes. It therefore it seems premature to conclude [^] that the data indicate a 
steeply-rising gluon in this x region. Artifically enhancing the gluon to fit the CDF low 
Pj' data, at the same time keeping the overall gluon momentum fraction fixed, inevitably 
depletes the large x part of the distribution, thus spoiling the fit to the WA70 data.! 

The hh cross section, being proportional to the square of the gluon distribution, is also of 
potential importance. Again, unfortunately, there are substantial theoretical uncertainties 
arising from a large scale dependence, and experimental uncertainties in reconstructing 
the total cross section from the measured distributions of B mesons, leptons, J/ip^s etc. 



in restricted kinematic ranges, see for example Ref. |42|. Also bb production at the FNAL 



pp collider samples x values {x < 0.01) just below the cross-over point where the various 
gluons have similar magnitudes, see Fig. 15. The most recent indications |^3|] are that 
while there is no serious disagreement between theory and experiment, the uncertainties 
are still much too large to allow any discrimination between gluon distributions at the 
level required by Fig. 15. 



In summary, there is still a considerable amount of uncertainty in the gluon distribution, 
particularly in the important medium-small and very small x regions. HERA will eventually 
provide information on the latter, and in the meantime large rapidity jet and prompt photon 
production at the FNAL pp collider are beginning to provide useful constraints. Current studies 
tend to use sets like MRS(Dq) and MRS(D'_) to quantify the effect of different small-x gluon 
shapes on cross section predictions, even though the corresponding quark distributions are now 
ruled out by the HERA F2 data. In a future study, we will explore the possibilities for varying 
the gluon distribution while maintaining a good fit to the structure function data. 



6. Hadroproduction of W bosons and top quarks 

In this section we present predictions for the total W and tt cross sections in pp collisions 
at 1.8 TeV. The former is being measured more and more precisely by both CDF and DO, and 
already provides a valuable cross check on the size of the u and d distributions. Interest in the 
latter has recently been rekindled by the evidence for top production presented by the CDF 
collaboration [J^, with the suggestion that the measured production cross section is larger than 



the theoretical prediction. 

(a) W production cross section 

While the W asymmetry probes the relative size of the u and d distributions, the total 
cross sections for W and Z boson production in pp collisions at ^/s = 1.8 TeV provide an 



^We are grateful to Steve Kuhlmann for discussions on this point. 
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important check of the overall magnitude of the quark distributions in a region of x where they 
are constrained by deep inelastic (mainly NMC and CCFR) data. Since the subprocess cross 
section is known to next-to-next-to- leading order [Q, there is very little theoretical uncertainty 
in the predictions once the parton distributions are specified. 

To illustrate this, we show in Fig. 19(a) the total W production cross section for the 
MRS(A,H,Dq,D'_) and CTEQ2M parton distributions. The central points are calculated as- 
suming Mw = 80.21 GeV/c^ and setting the renormalization and factorization scales equal to 
Mw- The errorbars show the effect of changing both scales simultaneously by a factor of 2. 
The first point to note is that the aw predictions for MRS (A), MRS(H) and CTEQ2M are very 
similar. (Notice also that the scale dependence is weak, and gives a similar variation in cross 
section as the spread in the three predictions.) This can be understood by recalling that aw 
is largely determined by the product u{xw, M^)d{xw, M^) where xw = Mw/y/s ~ 0.05. By 
comparing the relative size of the u and d distributions for the MRS (A), MRS(H) and CTEQ2M 
sets. Fig. 14, we can understand the similarity of the corresponding aw predictions shown in 
Fig. 19(a). First, we have that 

1iA)<m] ^ --(A) ^ .v,(H) . (37) 

The situation with CTEQ2M is more subtle. By comparing the u and d distributions at xw, 
we would conclude that 

d(A)>d(2M)} ^ -«'(A) > .„,(2M) . (38) 

whereas in fact the reverse is true. The reason lies in the different shapes of the W rapidity 
distributions. The CTEQ2M W cross section receives a significant contribution from 
produced at large positive (negative) rapidity, where a large-x w-quark interacts with a small-x 
d-quark. Especially for the former, the CTEQ2M distribution is significantly larger than the 
MRS (A) distribution, see Fig. 14, giving a distinctly different W rapidity distribution, as shown 
in Fig. 20 for pp + X. This is, of course, the origin of the much too large CTEQ2M 

lepton rapidity asymmetry discussed in section 4(b) and shown in Fig. 3. Comparing MRS(A) 
and MRS(H), we see that the distribution in yw for the latter is slightly larger for yw < 0, 
which is reflected in the total cross section shown in Fig. 19(a). Since the values and gluon 
distributions (for x ~ xw, see Fig. 18) for these three sets are very similar, the above orderings 
of u and d at Q = Mw simply reflect the orderings at the lower scales where the distributions 
are constrained by the fixed-target deep inelastic and asymmetry data discussed in the previous 
sections. Of the three sets only MRS (A) gives a consistent fit to all the data, and so we may 
conclude that aw (A) is the most accurate prediction. 

It is also worth commenting on the differences between the MRS(H), MRS(Dq) and MRS(D'_) 
W cross sections. The three parton sets have similar u- and d-quark distributions for x ~ xw 
and values characteristic of the fixed target deep inelastic scattering data [jlO|- There is, 
however, a rather large spread in the corresponding W cross sections shown in Fig. 19(a). These 
differences arise from the different gluon distributions in the three sets. We have seen in the 
previous section that x ~ 0.05 is precisely the region where the gluon is not well constrained. 
Thus even though the quarks are accurately determined at low Q for this x range, differences 



17 



in the gluon can affect the GLAP evolution to Q = and lead to differences in the quarks 
at this scale. This is graphically illustrated by the aw{H), aw(Do) and aw(D'-) predictions 
shown in Fig. 19(a). From Fig. 18, we see that the ordering of the gluons at low Q and x = 0.05 
is g(D'_) < 5'(H) < g(D'o), and hence the quarks of D'__/Dq increase more slowly /rapidly than 
those of H. This leads directly to (Jiy(D'_) < crwOi) < crvK(Do)- It would appear, therefore, that 
now that the quarks are very well constrained, it is the uncertainty in the gluon distribution 
around x ~ 0.05 which constitutes the largest theoretical uncertainty in the predictions for the 
W (and Z) cross sections at the FNAL pp collider. 

A precision measurement of aw at the percent level would obviously provide a useful addi- 
tional constraint on the distributions. The most recent values [ESl E6l for aw areP 



CDF(e,/i) : aw = 20.4 ± 0.4 (stat) ± 1.2 (sys) ± 1.4 (lum) nb 
DO(e, fi) : aw = 19.0 ± 0.6 (stat + sys) ± 2.3 (lum) nb (39) 

We should note that these cross sections relied on a luminosity measurement which has since 



been superceded by a new measurement of the total pp cross section by CDF If the 

new values are indeed some 10% higher, as suggested by the ratio of 'new' and 'old' total 
cross section measurements, then this would bring them more into line with the theoretical 
predictions shown in Fig. 19(a). However the errors on the data in Eq. (PD| ) are still much too 
big to provide the level of discrimination between the distributions discussed above. Indeed the 
accuracy of the theoretical prediction suggests that the W cross section could itself serve as a 
determination of the collider luminosity. 

Finally, we note that HERA will eventually be able to make a reasonably precise measure- 
ment of F2 in a region of Q intermediate between the fixed target data and Mw (see Fig. 17). In 
fact the first data in this kinematic region have recently been presented by HI p[ and ZEUS 0, 
although the errorbars are so large that no meaningful discrimination between the predictions 
is possible at present. 



(b) tt production cross section 



The evidence presented recently by the CDF collaboration |^ has focussed attention on the 
theoretical prediction for the top quark production cross section. Here we address the question 
of the uncertainty in the calculation of au coming from parton distributions. Fig. 19(b) shows 
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the next-to-leading order QCD predictions |^ for the total cross section for rrit = 174 GeV/c 
- the central value reported by CDF - with the factorization and renormalization scales set 
equal to m^. The errorbars on the predictions indicate the effect of changing the scales by a 
factor of 2. The relative spread in the predictions is somewhat less than for aw, since now the 
quark distributions are being probed in the range of x ~ 0.2 — 0.4 where the high-statistics 
BCDMS and CCFR data provide a very tight constraint.0 The more important point, however, 
is that unlike for aw, the dominant effect is the scale dependence. As discussed in detail in 



Ref. 1 50], this scale dependence is symptomatic of large contributions to the cross section from 



soft gluons at higher orders in perturbation theory. Fig. 21 shows the next-to-leading order 



^We have divided out the Standard Model leptonic branching ratio B{W — > Iv) — 0.108. 
^°Note that the qq annihilation process accounts for approximately 90% of the total cross section for this 
value of mt- 
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prediction for function of using the MRS (A) partons. The sohd curve corresponds 



to the scale choice Q 

m 



rrit and the dashed curves to Q = 2mt, hrrit. The data point is the CDF 



measurement 



(Ttt = 13.9 pb for mt = 174 ± 17 GeV/c\ 



7. Conclusions 

To determine the partonic structure of the proton we have performed a global next-to-leading 
order analysis of the data for deep inelastic scattering and related hard-scattering processes. 
We summarized the constraints from the various types of data in Table 1. We found that a 
surprisingly economical parametrization of the parton densities at = 4 GeV^ was able to 
give a satisfactory description of the data. We were able to include in the analysis for the first 
time significant new measurements of F2 at HERA, of the asymmetry A^y in Drell-Yan lepton 
pair production from pp and pn collisions, and of the rapidity asymmetry Aw at FNAL. 
The partons that we obtain - MRS (A) - are the only set which are consistent with all the data. 
(The only exception which remains is the incompatibility with the lowest x measurements by 
CCFR of the neutrino structure functions F2 and 0:^3.) 

The asymmetry data considerably tighten the constraints on the u, d, u and d distribu- 
tions. In particular, the NA51 measurement of A^y provides a missing link and pins down the 
combination d — u, which is essentially undetermined by the deep-inelastic structure function 
measurements. The measurements of the lepton rapidity asymmetry, A^/, provide a tight 
constraint on the ratio d/u with u and d having less influence. 

An interesting feature to emerge from the analysis is the link between the u and d distribu- 
tions at ~ and x ~ 0.05, which determine the W cross section at FNAL, and the u and 
d densities determined from the fixed-target data with ~ 20 GeV^. Even for x as large as 
0.05, the evolution from ~ 20 GeV^ to depends appreciably on the gluon distribution, 
particularly in the initial stages of the evolution. 

The measurements of F2 at HERA have opened up the small x regime, x ^ 10"'^. The 
dramatic growth of F2 observed with decreasing x, which had been anticipated in perturbative 
QCD models, is readily accommodated in the global fit. Indeed we found that the observed 
small X behaviour is entirely consistent with the precocious onset of BFKL leading log(l/x) 
dynamics, although alternative QCD explanations are not ruled out. 

We believe that our analysis has considerably improved the detailed knowledge of the quark 
densities. However, the same is not true for the gluon distribution. At present there are, with 
the possible exception of the WA70 prompt photon measurements, no reliable direct constraints 
on the gluon. In principle, bb and jet production cross sections depend directly on the gluon 
density, but ambiguities - due to scale dependence, jet recognition, higher-order effects, the 
choice of rrib etc. - mean that so far these data have provided little additional information on 
the gluon. Valuable indirect constraints on the gluon come from the momentum sum rule and 
from the observed structure of the deep-inelastic scaling violations. Our philosophy has been 
to keep the parametrization of the partons as economical as possible, and to only introduce an 
extra parameter when it is required to improve the description of the data. For this reason the 
parametrization of the gluon is remarkably simple, such that its behaviour is prescribed by the 
momentum sum rule and the optimum fit to the WA70 prompt photon data. It is therefore not 
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surprising that the gluons from the various sets of partons are quite similar above x ~ 10^^. 
Clearly this does not reflect the true ambiguity in the gluon distribution, but rather it is an 
artefact of our economical parametrization. We will attempt to address this problem in a future 
publication. We will also evolve our distributions to lower Q^, below Qq = 4 GeV^, and study 
the interplay of leading and higher twists, and hence provide a phenomenological representation 
of F2 at low Q^. 

For the future, we foresee a continued improvement in our knowledge of the partons. The 
increasingly precise measurements of F2 at HERA, particularly at low x and low [Q"^ ~ 
5 GeV^) may help to distinguish between BFKL and GLAP dynamics and may even reveal 
the onset of shadowing behaviour. So far the Drell-Yan asymmetry has only been measured 
at a: = 0.18. This leaves considerable ambiguity in the structure of d — u. The results of the 
asymmetry measurement by a forthcoming experiment at FNAL will cover an extended x 
interval and are eagerly awaited. The accuracy of the rapidity asymmetry measurements 
at FNAL will improve and thus provide extremely tight constraints on the behaviour oi d/u 
in the x region around 0.05. The neutrino measurements at low x are under further study by 
the CCFR collaboration and the discrepancy with the NMC measurements may be resolved. A 
change may effect the relative normalization between the data sets - an improvement here would 
be of great value. Finally, detailed studies of J/ip, cc, bb, prompt photon and jet production at 
FNAL and HERA, together with measurements of the longitudinal structure function F^, will 
help reduce the ambiguity in the gluon distribution. 
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Figure Captions 



[1] The HI and ZEUS data obtained from the 1992 HERA run, but pubhshed in 1993, 
together with earher data (interpolated to = 15 GeV^) obtained by the NMC |0 and 
BCDMS collaborations 0. The "pre-HERA" extrapolations obtained using MRS (D[),D'_) 
parton sets were expected to span the forthcoming HERA measurements. The "post- 
HERA" curve H is the result of a global analysis |jlO[ which included the HERA data. 



[2] The measurement of the asymmetry in Drell-Yan production in pp and pn collisions made 
by the NA51 collaboration at a; = y/r = 0.18. The curves obtained from MRS(H) 
and CTEQ2M ||lT| partons pre-date the measurement. The MRS(A) curve is obtained 
from the global fit, presented in this paper, which includes the NA51 data point. 

[3] The asymmetry A[y£) of the rapidity distributions of the charged leptons from W"^ — > £^z/ 
decays observed at FNAL as a function of the lepton rapidity y^. The curves are the 
next-to- leading order descriptions obtained using MRS(H) |jlO[, CTEQ2M |11| and the 



new MRS (A) partons. The MRS (A) analysis, presented in this paper, includes the data 
in the global fit. 

[4] The shaded band is the strange sea quark distribution, xs{x, Q"^ = A GeV^), determined by 
the CCFR collaboration from a next-to-leading order analysis of their dimuon produc- 
tion data. In addition they find that s/{u + d) is, to a good approximation, independent of 
X. Also shown is the MRS(A) input strange sea quark distribution, xs{x,Ql = 4 GeV^), 
with the 50% suppression, as given in Eq. (7). 

[5] The description of the EMC data ^Tj for F2 by the MRS (A) partons. We assume that 



the c — > /i -|- X branching ratio is 8%. 

[6] The description of the BCDMS and NMC measurements of the F2^{x, Q"^) structure 
function by the MRS (A) set of partons. The BCDMS data are shown with an overall 
renormalization by a factor of 0.98. 

[7] The continuous curves show the description of the NMC data for the structure func- 
tion ratio / given by the MRS (A) set of partons. The dotted curves are the 
predictions of the MRS(H) partons obtained from a global analysis which included an 
earher set of NMC data for F^/F^^. 

[8] The MRS (A) fit to the CCFR measurements p2| of the structure functions F2^ and 



xF^ . The data are shown after correction for the heavy target effects and after an 
overall renormalization of 0.95 required by the global analysis, see text. The statistical, 
systematic and heavy target correction errors have been combined in quadrature, see also 
Ref. i. 

[9] The MRS (A) partons shown as a function of x at Q'^ = 10 and 10^ GeV^. 

[10] The description of the preliminary ZEUS measurements of F2^ for x < 0.005 by the 
parton sets of Refs. |]^, |TD|. The long dashed curves show the fit to the data obtained by 
the MRS (A) analysis presented in the paper. 
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As for Fig. 10 but showing the prehminary HI measurements 0. 

x((J- m), as a function of x, at = 7 GeV^ obtained from the MRS(H) |0|, CTEQ2M 
|Tl| and the new MRS (A) sets of partons. 



A comparision at g2 = 20 GeV^ of the new MRS (A) partons of this analysis and the 
MRS(H) partons of Ref. pi. 



The u, d and d/u distributions at = m"^ obtained from the MRS(H) 0, CTEQ2M 
and the new MRS (A) sets of partons. 



The gluon distributions from sets Dq (dotted), D'_ (dashed) and A (continuous curve) at 
= 20 GeVl 

The gluon distributions and the corresponding predictions for dF2/d\ogQ'^ and Fl for 
the Dq, D'_ (dashed lines) and A (solid line) parton distributions, together with data from 
HI [0. 

The structure function at x = 0.05 as a function of Q, with data from NMC and 
the descriptions of the MRS(D'o, D'_, A, H) and CTEQ2M partons. 

The gluon distributions xg{x,Q^) of the MRS(D^,D'_,A,H) and CTEQ2M partons, at 
X = 0.05, function of Q. 

(a) Predictions for aw and (b) au (w^t = 174 GeV/c^) in pp collisions at = 1.8 TeV, 
from the MRS(Dq, D'„, A, H) and CTEQ2M partons. The central values are calculated 
using renormalization and factorizations scales equal to (a) Mw and (b) rrit, and the 
errorbars indicate the effect of increasing and decreasing these by a factor of 2. 

Next-to-leading order predictions for the rapidity distribution in pp collisions at 
y/s = 1.8 TeV from the MRS (A, H) and CTEQ2M partons. 

Next-to- leading order predictions for att at y/s = 1.8 TeV as a function of rrit using 
MRS (A) partons. The curve is calculated using renormalization factorization scales equal 
to rrit, and the band corresponds to changing scales to ^m^ and 2mt. The data point is 
from the CDF collaboration . 
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